Spurious fixed points in frustrated magnets 
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We analyze the validity of perturbative estimations obtained at fixed dimensions in the study of 
frustrated magnets. To this end we consider the five-loop /3-functions obtained within the minimal 
subtraction scheme and exploited without e-expansion both for frustrated magnets and for the well- 
controlled ferromagnetic systems with a cubic anisotropy. Comparing the two cases it appears that 
the fixed point supposed to control the second order phase transition of frustrated magnets is very 
likely an unphysical one. This is supported by the non-Gaussian character of this fixed point at 
the upper critical dimension d = 4. Our work confirms the weak first order nature of the phase 
transition and constitutes a step towards a unified picture of existing theoretical approaches to 
frustrated magnets. 
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Introduction. Although undoubtedly successful to de- 
scribe the critical behavior of 0(7V)-like models, per- 
turbative field theory is still unable to provide a clear, 
non-controversial understanding of the physics of more 
complex models among which are the famous Heisen- 
berg or XY frustrated magnets (see and references 
therein). At the core of the problem is that different 
kinds of perturbative approaches, performed up to five- 
or six-loop order, lead to contradictory results: first order 
phase transitions are predicted within the e (or pseudo- 
e)- expansion |^ Q whereas a second order transi- 
tion is found in the fixed-dimension (FD) perturbative 
approaches performed either in the minimal-subtraction 
(MS) scheme without e-expansion ^ or in the massive 
scheme . In fact FD results for frustrated magnets are 
neither supported by experiments nor by Monte Carlo 
simulations S S ^1 (see however p where a scal- 
ing behavior is found). They also disagree with the pre- 
dictions obtained from the non-perturbative renormaliza- 
tion group (NPRG) approach [l|,|ll| that predicts (weak) 
first order phase transitions in d = 3 in agreement with 
the e-expansion analysis. 

In this letter we shed light on the discrepancies encoun- 
tered in perturbative approaches to frustrated magnets 
by showing that the FD approaches very likely lead to 
incorrect predictions as for the critical physics in three 
dimensions. Our key-point is easy to grasp already for 
the simplest — 0{N) — model. In this case, the (non- 
resummed) renormalization group (RG) /3-function at L 
loops is a polynomial in its coupling constant u of order 
L + 1. Thus, it admits L + 1 roots u*, f3{u*) = 0, that 
are either real or complex. Within the e-expansion ap- 
proach, the only non-trivial fixed point (FP) retained is 
such that w* ~ e where e = A — d. On the contrary, in the 
FD approaches, no real root can be a priori discarded. 
As a result, the generic situation is that the number of 
FPs as well as their stability vary with the order L: at a 



given order, there can exist several real and stable FPs 
or none instead of a single one. This artefact of the FD 
approach is already known and was first noticed in the 
massive scheme in d = 3 The way to cope with it 

is also known: resumming the perturbative expansion of 
/3(u) (see e.g. 13, 14]) is supposed both to restore the 
nontrivial Wilson-Fisher FP and to suppress the non- 
physical or "spurious" roots. This is indeed what occurs 
for the 0{N) model for which the FP analysis performed 
on the resummed /3-function of FD approaches enables to 
discriminate between physical and "spurious" FPs. We 
argue that the situation is very different for frustrated 
magnets. Indeed, considering the /3-functions derived at 
five loops in the MS scheme and using the same resum- 
mation procedure as Calabrese et al. 0| we show that 
the FP found in d = 3 without expanding in e is spurious 
although it persists after resummation. Our conclusion 
is based on two main facts: (i) analyzed with the same 
FD approach a similar FP is found for the ferromagnetic 
system with cubic anisotropy in contradiction with its 
well established critical physics (ii) there are strong in- 
dications that these FPs survive in the upper critical di- 
mension d = A where they are found to be non-Gaussian 
which raises serious doubts on their actual existence. 

Resummation. To investigate the five-loop f3 functions 
derived in the MS scheme we resum them, following 
using a conformal mapping Borel transform suitable for 
series involving two coupling constants ui and U2. We 
consider any function / of ui and U2 as a function of ui 
and z = M2/U1, supposed to be known through its series 
expansion: 



f{ui,z) =^fn{z) 



(1) 



An important hypothesis underlying the procedure of re- 
summation used here and in Q is that one can safely 
resum in the ui direction while keeping z fixed, i.e. 
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without resumming with respect to U2. Under this hy- 
pothesis a resummed expression associated with / reads: 



(a, a(z), 6; z) j dt 





[uj{uit-z)f 
[1 - uj{uit; 2)]°" 

(2) 

with uj{u] z) — (^1 + a{z) u — + a{z) u + 1) and 

where the coefficients dn{a,a{z),b] z) are computed so 
that the re-expansion of the r.h.s. of in powers of 
Ml coincides with that of J^l- In principle the parame- 
ters a(z), b and a are determined (i) by the asymptotic 
behavior of /„: /„_»oo ~ (— a(z))"n!n^ and (ii) by the 
strong coupling behavior of /: /(wi — > oo,z) ~ u"^^. 
However in general only a is known and the other pa- 
rameters are considered as free or variational. 

Frustrated magnets. The Hamiltonian relevant for 
frustrated systems is given by: 



U2_ 
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(3) 

where the 4>i are iV-component vector fields. We apply 
the resummation procedure described above without e- 
expansion J^] to the I3ui functions, i = 1, 2, obtained at 
five loops in the MS scheme |^. More precisely, as in 



3ui{ui^z) + eui)/ui , i = 1,2, instead 



0, we resum 

of Pui{ui,z) [26j. In this model the region of Borel- 
summability is given byui — U2/2>0 while a{z) = 1/2 
and b and a are typically varied in the ranges [6, 16] and 
[—0.5, 2]. As in all other approaches, one finds, in agree- 
ment with , that there exists a curve (parametrized by 
Nc{d) or its reciprocal dc{N)) such that for d < dc{N) a 
stable FP C+ governs the critical properties of the sys- 
tem. The curves Nc{d) obtained within this scheme, 
the NPRG approach and the e-expansion are given in 
Figffl This figure sums up the discrepancies between 
these approaches. The S-like shape of the curve N^^{d) 
obtained within the perturbative FD approach is such 
that C_|^^ exists for all > 2 in d = 3 contrary to 
the other approaches in which a FP C+ exists only for 
AT > Afe(d = 3) ~ 5 Q. 

The main question is to know whether one can trust 
the results obtained at FD, in particular the existence of 
FPs in d = 3 for A^ = 2,3. To answer this question we 
examine the convergence of the FD perturbative results 
through their sensitivity with respect to variations of the 
resummation parameters b and a as well as the order L 
of the computation. We concentrate our study on the 
(real part of the) stability critical exponent uj at the FP 
which is focus 0]. We have checked that similar 
behaviors are obtained for the exponent u 0]. In prac- 
tice, following we optimize uj{a,b,L) by choosing a 
such that Lu{L -I- 1) — uj{L) is minimal and b such that 
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FIG. 1: Curves Nc{d) obtained within the e-expansion (A'^^), 
the MS scheme without e-expansion [N^^) and the NPRG 
approach (A^^^^''). The resummation parameters for the MS 
curve are a = 1/2, 6 = 10 and a = 1. The part of the 
curve N™ below S corresponds to a regime of non-Borel- 
summability. 



tj(6) is (quasi-)stationary. In FigEl we display uj{b) in 
the Heisenberg case for L = 4 and 5 for typical values of 
a. At a given order the variations of u with b are small 
around the stationary point and do not exceed 30% from 
four to five loops in agreement with [^. Similar results 
are obtained in the XY case. From this analysis alone, 
one could conclude that a FP indeed exists in d = 3 al- 
though with error bars on the critical exponents far larger 
than those obtained for the Ising model with the same 
methodology (much less than 1% at this order). Let us 
now study the cubic model along the same lines. 
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FIG. 2: The (real part of the) critical exponent as a function 
of b at five (upper curves) and four (lower curves) loops for 
a = —0.5, and 0.5 for the frustrated model (A'^ = 3). 

The cubic model. We now consider the ferromagnetic 
model with cubic anisotropy whose Hamiltonian is: 



n 



2J.2 



m 4> 



u 
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'4>\ 



N 

i=i 



(4) 

with <j> a TV-component vector field. The Hamiltonian Q 
is used to study the critical behaviors of numerous mag- 
netic and ferroelectric systems with appropriate order pa- 
rameter symmetry (see e.g. 0,0]). The /3- functions are 
known at five-loop order in the MS scheme 'io] and at 
six-loop order in the massive scheme Apart from 
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the Gaussian and an Ising FP, there exist two FPs: the 
0{N) symmetric FP {u* ^ 0,v* = 0) and the mixed one 
M {u* ^ 0,v* ^ 0). The Ising and Gaussian FPs are 
both unstable for all values of N. The 0{N) FP is sta- 
ble and M is unstable with w* < for < Nc and the 
opposite for N > Nc- Nc has been found to be slightly 
less than 3 [Ullli. 

Let us now analyze the FP structure of the model 
within the MS scheme without e-expansion by ap- 
plying the conformal mapping Borel transform The 
parameter a(z = v/u) entering in ^ is now given by 
a{z) = 1 + z for 2; > and a{z) = 1 + z/N for z < while 
the region of Borel-summability is given by the condition 
u+v > and Nu+v > 0. Within this scheme one surpris- 
ingly observes that, in addition to the above mentioned 
usual FPs, there exist, in a whole domain of parameters 
b and a, several other FPs that have no counterpart in 
the e-expansion. In particular, one of them that we call 
P (which is stable and such that u* > 0,v* < 0) exists 
for any value of < 7.5 and lies in the region of Borel- 
summability u + v > 0. The presence of this FP, if taken 
seriously, would have important physical consequences 
since it would correspond to a second order phase transi- 
tion with a new universality class. However as far as we 
know, no such transition has ever been reported and the 
existence of P has to be considered as an artefact of the 
FD analysis. Note finally that P is found to be a focus 
FP, a striking similarity with the frustrated case. 

At this stage, it is very instructive to perform for P 
the same convergence analysis as we have done for C+ 
in the FD analysis. In Fig^we give uj{b) for i = 4, 5 
and for three different values of a. Surprisingly, we find 
results that display similar convergence properties as in 
the frustrated case (with a spreading of the results of 
order 40%). 

This study shows that contrary to what one could 
naively deduce, the convergence properties displayed in 
FigEldo not characterize P as a genuine FP. Comparing 
the studies of the frustrated and cubic cases one is led to 
the same conclusion for Cj^ for A^ = 2,3ind = 3. Thus, 
in order to conclude whether C^^ is a genuine FP or an 
artefact of the FD analysis, one has to resort to another 
criterion. 

Fixed point at the upper critical dimension. To do 
this let us perform the following observation: the choice 
of (f>^-\ike Hamiltonian in Eq. Q relies on the physical hy- 
pothesis that the upper critical dimension is four. From 
the perturbative point of view, this means in particu- 
lar that the theory is infrared trivial — controlled by 
the Gaussian FP — in four dimensions. Thus, by con- 
sistency, one must retain among all solutions of the FP 
equations in d = 4 — e only those that are at a distance 
of order e of the Gaussian FP. Therefore a practical way 
to check whether a FP found at a given dimension, d = 3 
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FIG. 3: The (real part of the) critical exponent a; as a function 
of b at five (upper curves) and four (lower curves) loops for 
a = 1, 1.5 and 1.7 for the cubic model (N=2). 

for instance, is a genuine FP or is just an artefact of per- 
turbation theory is to follow it by continuity up to the 
upper critical dimension 0, . If the FP survives as a 
non-Gaussian FP at this dimension this is a signal that 
it is spurious. 

Let us apply this criterion simultaneously to P and 
. We present our results in Fig2] where we have 
displayed the coordinates u* and u* associated to these 
FPs as functions of d. Manifestly, they both survive ev- 
erywhere above d = 3 and are not Gaussian in c? = 4. In 
the cubic case this is true for all values of N for which 
P exists while, in the frustrated case, this happens for 
all values of iV < 6. Thus, according to our criterion, P 
is, as expected, always found to be spurious. As for 
for N < 6 and, in particular, for TV = 2, 3 our criterion 
strongly suggests that it is also spurious in d = 3 27*1. 
Note that the FPs P and Cj^ lie out of the region of 
Borel-summability in d = 4. Thus their coordinates can- 
not be determined accurately although their existence as 
non-Gaussian FPs is doubtless. 

u* , Ui* 
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FIG. 4: The u* coordinate of the FP P (TV = 2, upper curve) 
and the ul coordinate of the FP C™ (N — 3, lower curve) 
as functions of d. 

Conclusion. It appears from our study that the FP 
identified in the FD approach is very likely spuri- 
ous. The transition should thus be of — possibly weak — 
first order in agreement with NPRG and e-expansion ap- 
proaches. It remains to explain the failure in the resum- 
mation procedure used in the FD approach, Eq.|21). As 
already emphasized this procedure relies on the hypoth- 
esis that resumming with respect to mi, keeping a poly- 
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FIG. 5: Four curves Nc^{d) for different values of the param- 
eter a and the curve N^{d). From right to left a — 0.5, 1.3, 1.5. 
The parameters b and a are 6 = 10 and a — 1. The parts 
of the curves below the black dots correspond to a regime of 
non-Borel summability. 



nomial structure in U2, is sufficient. Alternatively a re- 
sumniation of the series with respect to the two coupling 
constants could be required to obtain reliable results (see 
for instance '23^1 for the randomly diluted Ising model). 
Postponing these considerations for a future publication 
jT3| we assume that the use of Eq.jSJ as such is justi- 
fied. Then a possible origin of the failure in the resom- 
mation procedure could be that the series considered are 
not known at large enough order to reach the asymptotic 
behavior. In this case there would be no reason to fix the 
parameter a at its asymptotic 1/2 value and one should 
have to vary it as b and a to optimize the results [l8j |. 
We display in FigElthe curves N^^{d) for different val- 
ues of a. The part corresponding to large values of Nc, 
typically for Nc ^6, is almost insensitive to the varia- 
tions of a whereas this is clearly not the case for smaller 
values of Nc- In particular, for sufficiently large values 
of a, typically a > 1.5, the S-like part is pushed below 
d = 3 so that N^{d), Nj^{d) and N^^^^{d) are then 
compatible everywhere for 3 < d < 4. Thus, under our 
hypothesis, all qualitative differences between the differ- 
ent approaches disappear so that the problem would 
boil down to a question of order of computation. It is 
clear that our heuristic argument, although appealing, 
requires more rigourous developments |l7j . 

Note finally that our present considerations surely per- 
tains to the case of frustrated magnets in d = 2 
Indeed we have checked that the FP found in d = 2 is 
continuously related to the FP in d = 3 which makes 
its existence doubtful. Our conclusions could also apply 
in other situations where FPs that have no counterpart 
in the e-expansion is found, as it is the case, for instance, 
in QCD at finite temperature [j^ . 

We wish to thank P. Azaria, P. Calabrese, R. Guida 
and J. Zinn- Justin for useful discussions. 
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